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MAGNETOGASDYIMAMIC COMPRESSION OF A COAXIAL PLASMA ACCELERATOR 
FLOW FOR MICROMETEOROID SIMULATION 


I. INTRODUCTION 


The extended duration of manned space flight missions and unmanned operations 
with very sensitive equipment makes it necessary to provide an adequate protection 
against micrometeoroid impact Experimental work on the meteoroid impact problem is 
performed using several types of projectile accelerators which provide high velocity 
simulated micrometeoroid fluxes in the laboratory 

The light gas gun is the most commonly used tool. It accelerates metallic and 
nonmetalhc particles with diameters from 0 015 cm to 0.635 cm to velocities as high as 
9.5 km/s, and the maximum particle energy is about 1 4x 10^ J (0 316 cm diameter 
lexan cylinder at 9.5 km/s) [1,2]. More recently, detonation waves in explosives have 
been used to generate a high speed gas flow [3] which is used to accelerate a projectile 
For small particle diameters, 1 micron and below, electrostatic accelerators are very 
useful and reliable micrometeoroid flux simulators [4,5] The energy of the accelerated 
particles js limited by the charge-to-mass ratio Very small particles with diameters of 
about 0 1 micron can be accelerated to velocities of approximately 70 km/s. Several 
efforts have been made to transfer the energy stored m a capacitor bank to particles for 
micrometeoroid simulation The most successful accelerator has been developed by Scully 
[6] His device can accelerate particles with diameters of approximately 50 microns to 
velocities of about 20 km/s. 

The expenmental and theoretical investigations described in the subsequent 
sections of this report were aimed at the development of a new concept for the 
generation of a high density plasma for particle acceleration The energy stored in a 
capacitor bank is used to generate a high velocity, hot plasma in a coaxial accelerator 
[7,8,9,10,11] The plasma leaving the coaxial accelerator is compressed in a compressor 
coil which IS mounted on the exit end of the accelerator The most interesting feature of 
this new combination of coaxial accelerator and compressor coil is that the compressor 
coil does not need an external power source It obtains the energy needed to compress 
the plasma from the capacitor bank which is discharged across the coaxial accelerator 
The new device is, therefore, very simple and can be added easily to existing coaxial 
plasma accelerators It has been demonstrated that glass beads of 125 to 700 micron 
diameter can be accelerated to velocities of 1 1 km/s to 5 km/s, respectively The 
accelerator configuration is presently being investigated to determine the optimum 
performance cntena 



II. DESCRIPTION OF THE PLASMA HYPERVELOCITY ACCELERATOR 


The plasma hypervelocity range for micrometeoroid simulation at the George C 
Marshall Space Flight Center’s Space Sciences Laboratory consists of a capacitor bank 
with switch assembly, the test range with plasma accelerator, and installations for plasma 
diagnostics and particle velocity measurement 


A. Capacitor Bank and Switch Assembly 

The capacitor bank consists of 72 individual 1 5-microfarad capacitors connected 
m parallel in groups of 18 capacitors All of the expenments were made with only one 
quarter bank The maximum charge voltage was 20 kV and the corresponding stored 
energy was 54 kJ. Three capacitors connected m parallel by copper bus bars form a 
module, and every module is connected to the switch assembly by two coaxial cables All 
the cables are the same length. Each quarter bank is equipped with four ignitron switches 
in parallel. These are fired by a trigger assembly containing a set of capacitors and a 
thyratron to prevent back current to the igniters of the ignitrons. The switch assembly is 
connected to the coaxial accelerator by six coaxial cables. 

The capacitor bank is charged from a high voltage supply which is controlled by 
an auto transformer and is disconnected before firing For emergencies a dump circuit is 
incorporated to discharge any energy stored in Ihe capacitor bank through the charging 
network within 30 seconds at a peak current of 4 microamps If a prefire occurs, the 
relays connecting the high voltage power supply and the capacitor bank become 
de-energized The trigger assembly has its own high voltage power supply and is charged 
to 5 kV It IS activated by a 150-V pulse from a delay generator which, in turn, is 
triggered either by an external input or manually A schematic diagram of the capacitor 
bank, charging circuit, switch assembly, and test range is shown m Figure 1 Table 1 lists 
the characteristics of the equipment 


B. Test Range with Plasma Accelerator 

The test range is made of a 2 0-m long, 1 5-cm diameter aluminum tube with a 
flange at both ends (Fig 2) It has one 38- by 12 5-cm window of 2 5-cm thick Plexiglas 
for photographing the plasma flow and the particle acceleration The inner surface of the 
window IS covered with a 3-mm replaceable Plexiglas cover for protection against the hot 
plasma. Two other windows on top of the range have 3. 5-cm thick aluminum covers 
with 5 0-cm diameter, vacuum-tight, feed-through plugs A vacuum pump is connected to 
the range tnrough a flexible hose The coaxial plasma accelerator with compressor coil is 
mounted m a fiberglass plate which closes the front end of the test range This plate has 
two vacuum-tight coaxial cable feed-through plugs for the measurement of the 
compressor coil current The other end of the test range is closed by an aluminum cover 
with another feed-through plug and a vacuum gage. 
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Figure 1 Schematic diagram of plasma hypervelocity range 


The coaxial accelerator with compressor coil is shown in Figure 3. The center 
electrode (1), the insulating plug (2), and the annular electrode (3), form the coaxial 
accelerator. An aluminum foil (4) is placed between the center electrode and the annular 
electrode before every finng. The positive side of the charged capacitor bank is connected 
to the center electrode through the switch assembly, and the negative side is connected 
to the annular electrode An insulating nng (5) separates the compressor coil (6) from the 
annular electrode and also serves as a mount The front (narrow) end of the compressor 
coil IS connected to the potential of the annular electrode by the holder rail (7) and is 
held in position by the coil spacer plate (8) The coil spacer adapter (9) keeps the holder 
rails in position In the increasing inductance mode configuration, the insulation ring (5) 
was not used and the holder rail (7) was disconnected from the annular electrode 


C. Diagnostics and Particle Velocity Measurement 

Photographs of the plasma flow and the accelerated particles were made through 
the 38- by 12.5-cm window in the range Two single frame image converter cameras were 
used for high speed photographs They were triggered by the pulse from the delay 
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generator to the capacitor bank firing circuit Using the adjustable exposure delay of the 
cameras, the time interval between the two exposures could be set between 0 and 100 
microseconds 


TABLE 1 PHYSICAL CHARACTERISTICS OF EQUIPMENT IN COAXIAL 

ACCELERATOR CIRCUIT 


— 

Equipment/Characteristic 

Description 

Capacitors 


Peak dc Voltage Rating 

20 kilovolt 

Capacitance 

15 microfarad +20%-! 0% 

Inductance 

40 nanohenry max 

Energy 

300 joule 

Voltage Reversal 

90% 

Peak Current 

2x10^ ampere (direct short) 

Life Expectancy 

10^ discharges 

Cables 


Type 

20 P 2 

Peak Voltage 

20 kilovolt 

Reversal Level for Full Nominal Life 

19 kilovolt 

Nominal Life 

10® pulses 

Characteristic Impedance 

16 ohms ± 1 ohm 

Inductance (at 200 kHz) 

105 ± 6 6 nanohenry /meter 

Delay 

6 56 nanosecond/meter 

Normal dc Resistance 


Inner Conductor 

0 568 ohm/kilometer 

Outer Conductor 

3 335 ohm/kilometer 

Minimum Bending Radius 

12 7 centimeter 

Capacitance 

407 picofarad/meter 

Ignitrons 


Type 

GL-7703 

Peak Anode Voltage (Forward) 

20 kilovolt 

Peak Anode Voltage (Reverse) 

20 kilovolt 

Peak Anode Current for One-Half-Cycle of 20 microsecond 

10® ampere 

Tube Inductance 

4 X 10'® henry 

Ionization Time 

0 5 microsecond 

Ignitor Ratings 


Ignitor Voltage 


Forward Open Circuit 

1500 volts min 


3000 volts max 

Inverse Voltage 

5 volts max 

Peak Ignitor Current 

250 ampere max 


200 ampere min 

Length of Firing Pulse, Sine Wave 

5 microsecond mm 


10 microsecond max 
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PHOTOMULTIPLIER ASSEMBLY 



Figure 2 Test range with plasma accelerator. 



Figure 3. Coaxial accelerator with compressor coil 
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An eight-camera unit made up of eight individual image converter cameras was 
used for sequential photographs In this unit the plasma image is deflected by a pnsm, 
which IS located m the center of a drum-shaped container, to the eight image converter 
tubes located at 45-deg intervals. For particle detection m the very bright plasma light, 
filters were placed between the prism in the center and the image converter tubes. The 
eight cameras could be triggered at a maximum time interval of 1 microsecond between 
successive frame exposures, whereas, the minimum exposure time was 10 nanoseconds 
and the maximum exposure time was 100 nanoseconds The trigger pulse for the image 
converter cameras came from a transient light detector or from the delay generator, 
which also tnggered the capacitor bank firing circuit In both cases, the trigger pulse for 
the first camera was recorded on an oscilloscope 

Light probes were placed m vanous positions m front of the plasma range window 
to measure the plasma front propagation velocity The output of these detectors was 
recorded on an oscilloscope The particle velocity was measured with the eight-camera 
unit and an impact detector. Using light filters, the particles could be photographed 
when they left their position at the Mylar foil. The time delays between the eight 
photographs were recorded on time interval counters. These delay times, together with 
the recorded trigger pulse for the first camera, bracketed the time when the particles had 
left their position on the front of the compressor coil 

The impact detector (Fig. 2) placed between 0 2 and 0.4 meters downrange had a 
detector mounted coaxially with the plasma accelerator and compressor coil An 
aluminum sheet was held m place by a nut The accelerated particles perforated this 
sheet, and a light flash was generated which was deflected by a mirror into a 
photomultiplier mounted on the top window cover of the range A hole was punched 
into one aluminum sheet to determine the origin of the light recorded by the 
photomultiplier. The plasma light coming through this hole for a discharge with no beads 
had the same intensity as the signal from a perforation by an accelerated particle It may, 
therefore, be concluded that the recorded light flash from an aluminum sheet perforation 
by an accelerated particle is generated by the plasma, either m the accelerator or m a 
shock wave m front of the impact detector 

The output of the photomultiplier detector was recorded on the same dual beam 
oscilloscope as the trigger pulse for the first camera The particle flight time, together 
with the known distance between the compressor coil front end where the photographs 
were made and the aluminum sheet in the particle impact detector, was used to 
determine the average particle velocity The capacitor bank voltage was measured with an 
inductive pickup cod mounted on top of the positive leads of the coaxial cables leading 
from the ignitron switch assembly to the coaxial accelerator The shields of the coaxial 
cables were stripped from the cables in the switch assembly, the positive leads go through 
the ignitrons and the shield goes through a ground plate The compressor coil current was 
measured indirectly by a pickup coil m front of a loop m the backstrap lead from the 
compressor coil to the center electrode To eliminate electromagnetic disturbances, a 
second pickup coil was mounted on the other side of the coaxial accelerator to record 
only the electromagnetic disturbances Both signals were then subtracted m a differential 
amplifier and then recorded on an oscilloscope 
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III. PLASMA COMPRESSION AND PARTICLE ACCELERATION 


Glass beads with diameters between 0.125 mm and 0.795 mm were used for 
acceleration expenments There was no indication of any forces other than gasdynamic 
forces acting upon these beads 


A. Drag Acceleration of Particles for Micrometeoroid Simulation 

In gasdynamics, the force acting upon a body may be expressed as the sum of a 
dynamic pressure term and the product of the pressure gradient times a characteristic 
length AZ, both multiplied with a charactenstic surface S, which is assumed to be the 
same here for both terms 


F = (Cd y V" + VPAZ) S 


( 1 ) 


For constant particle size and flow velocity, V, the dynamic pressure term is only a 
function of the flow density p Plasma accelerators for neutral particle acceleration were, 
therefore, designed for high flow densities 

In a magnetogasdynamic flow the pressure gradient term VP may become large if 
the plasma density is high enough to permit an immediate transfer of the electromagnetic 
force into gasdynamic pressure Maximum gasdynamic acceleration of neutral particles, 
therefore, is obtained by high density flow and a large negative pressure gradient. The 
density of the plasma leaving the coaxial accelerator was measured with the drum probe 
shown in Figure 4 A metallic foil is attached to a Plexiglas housing The probe in front 
of the foil indicates the arrival of the electrically conducting plasma, and a current flows 
between the potential drop across the resistance R1 in Figure 4 The nse-time of the 
circuit shown m Figure 4 is 1 nanosec/Volt, with the capacitance C in parallel with the 
battery B. The plasma pushes the foil towards the center electrode. Upon contact, the 
current through R1 is cut short and the potential drop across R1 becomes zero The time 
dunng which a potential can be measured across the resistance R1 is the time during 
which the foil moves from the initial position to the center electrode. 

From the acceleration of the foil with the mass M and the law of the 
conservation of momentum of the plasma flow hitting the foil, the density p of the 
plasma flow can be calculated Using the data from Table 2, the following value is 
calculated [7] 


p = 2 5 X 10'^ kg/m^ 


( 2 ) 
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FOILMOVING FOIL HITTING 

CENTER ELECTRODE 



DRUM PROBE 

Figure 4 Drum probe for impulse measurement. 
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TABLE 2. MEASURED DATA FOR TEST FACILITY 


Capacitor Bank Charging Voltage 

Uq = 16 kV 

Initial Pressure in Test Range 

Pq = 5 X 10'^ Ton- 

Foil Mass 

DC 

O 

II 

Displacement of Foil 

S = 0.35 cm 

Foil Radius 

R = 1 25 cm 

Acceleration Time 

T = 6.0 X 10'« s 


Another estimate can be obtained by dividing the mass of the alum mum foil inserted into 
the coaxial accelerator by the volume of the coaxial accelerator 


^FOIL 

''^COAX 


1 06 X 10'* kg/m^ 


(3) 


A comparison of the results from equations (2) and (3) shows that the measured density 
is approximately one-fourth of the calculated value This indicates that only a part of the 
aluminum foil mass is accelerated High-speed photographs using a movie camera show 
that plasma columns are accelerated at every half-period of the circuit consisting of the 
capacitor bank, cables, and the coaxial accelerator 

The measured density can be introduced into equation (1) Assuming that there is 
no contnbution from the pressure gradient term, the acceleration of a glass bead may be 
computed using the dynamic term Fp of equation (1), 


a 




/ ^PLASMA 
\ ^GLASS 


.4 R 




BEAD, 


(4) 


With the drag coefficient Cq = 20 from the Newtonian approximation for hypersonic 
flows and the glass density p = 2 5 g/cm^. 
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a 


3 V2 

Tr 


X 10"^ m/sec^ 


(5) 


For V = 40 km/s and R = 0.5 mm, the acceleration is 


a = 24 X 10’^ m/sec^ 


( 6 ) 


Using the measured acceleration time T of the foil for the projectile acceleration, the 
final velocity of the accelerated glass bead becomes 


V = a T = 144 m/sec 


(7) 


This result was venfied by the experiment No measurable acceleration was obtained 
when the beads were attached to a Mylar foil and placed across the end of the coaxial 
accelerator Therefore, an attempt was made to find a region of higher density inside the 
coaxial accelerator. 


B. Coaxial Accelerator for Particle Acceleration 

The electromagnetic force in the coaxial accelerator without a compressor coil is 
generated by the interaction of the radial current density jj. between the center electrode 

and the annular electrode, and the azimuthal magnetic induction around the center 
electrode Both are functions of the radial distance r from the center electrode, 

Jr ~ 7 and ~ 7 (8) 

The axial electromagnetic force is 


■’r ~ 


(9) 
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and the axial momentum equation becomes 


d V 


d t 


z 3 P 

+ = F. 


3 z 


( 10 ) 


The maximum acceleration and pressure gradient will occur m the vicinity of the 
center electrode. The plasma leaves a coaxial accelerator with equally long center and 
annular electrodes in a conical shape, as shown in Figure 5. Dunng the acceleration, 
plasma is forced radially toward the wall of the annular electrode. After 10 to 20 firings, 
the inner surface of the annular electrode is covered with aluminum deposits. The 
maximum acceleration of the particles was observed when they were placed near the 
center electrode, as shown in Figure 6 A copper, nng-shaped holder was inserted into the 
annular electrode, and the center electrode was shortened If the nng holder is in- metallic 
contact with the annular electrode, then the secondary discharge indicated m Figure 7 
occurs. This discharge was photographed with the center electrode protruding out of the 
annular electrode (Fig. 8). The discharge starts when the plasma fills the gap between 
the center electrode and the ring holder The azimuthal magnetic field around the 
discharge mteracts with the axial current density and produces a radial force and an 
axially directed force due to the conical shape of the discharge A thin Mylar foil with 
glass beads attached to the foil by a thin film of vacuum grease was placed across the 
nng holder The maximum particle velocity (1.0 to 2.0 km/s with 0.7 mm diameter 
beads) was achieved with the holder in front of the center electrode inside the annular 
electrode. 


INSULATOR 




PHOTOGRAPH 
(5 ns EXPOSURE TIME) 


Figure 5 Plasma leaving the coaxial accelerator. 
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HOLDER 
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PLASMA 

FLOW 


SECONDARY 

DISCHARGE 


Figure 6. Particle acceleration near center electrode of coaxial accelerator. 



Figure 7. Description of secondary discharge between center electrode 

and copper ring 


It should be noted that the configuration with the center electrode longer than 
the annular electrode (Fig 7) produced a very high velocity luminous front travelling 
down the range if no Mylar foil was placed on the ring holder The plasma arrived at the 
gap between the center electrode and the ring holder with a velocity corresponding to the 
coaxial accelerator configuration. The secondary. discharge provides additional acceleration 
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RING 


Figure 8 Photograph of secondary discharge 

because of its conical geometry and, thus, acts as a second acceleration stage Finally, the 
coaxial accelerator without the compressor coil produces a fast, low-density plasma The 
vaponzed and ionized aluminum foil does not stay together, and the plasma leaves the 
accelerator as a column which is at least as long as the accelerator 


C. Compressor Coil Configuration 

It IS more effective to use a compressor coil, as described earlier and shown m 
Figure 3, mounted on the end of the coaxial accelerator. The emerging plasma initiates a 
discharge between the center electrode and the copper compressor coil, which is 
connected to the annular electrode. If the spacing between the coil windings is large 
enough, then the electnc current will flow through the copper coil windings and wdl 
generate inside the coil a mostly axial magnetic field (Fig. 9). This magnetic field changes 
with time and induces an azimuthal current density in the plasma inside the compressor 
coil which interacts with the magnetic field and results m a radial force on the plasma 
While the axial magnetic field is increasing, this radial force acts toward the axis, keeping 
the plasma inside the compressor coil (Fig. 10). The current from the center electrode 
to the compressor coil has axial and radial components. The azimuthal magnetic field and 
the radial current density component generate an axial force as m the coaxial accelerator 
Therefore, the plasma will be compressed and accelerated into the compressor coil at the 
same time. The current density from the center electrode to the compressor coil is 
controlled by a local equihbnum between the radial force, which tends to push the 
plasma away from the compressor coil surface and the current density needed to produce 
this force 
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INDUCED CURRENT IN PLASMA. 


Figure 9. Compressor coil with magnetic field and plasma current density 


CURRENT FROM COMPRESSOR COIL 



Figure 10 Photograph of plasma flow m compressor coil during 

compression period. 
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There are two ways to connect the compressor coil to the annular electrode 
potential and they lead to two different modes of operation the decreasing and 
increasing inductance modes. In the decreasing inductance mode, the compressor coil is 
insulated from the annular electrode by an insulator nng (Fig 3) The end of the 
compressor coil is connected to the annular electrode through the holder rail As the 
plasma moves into the compressor coil, the current density will assume some distribution 
along the compressor coil The current flow to the compressor coil wide end, which is 
close to the coaxial accelerator end and attached to the insulator ring, will have to 
traverse all the windings of the compressor coil, whereas the current going to the narrow 
end of the compressor coil will flow only through the final winding In an equivalent 
circuit for this configuration every current element must be combined with a different 
inductance, varying from the total inductance L] at the wide end to a small inductance 

at the narrow end. If AC is a length element of the compressor coil, then the position 
measured from the wide end of the compressor coil is given by lAC, with i varying from 
zero to n, where nAC = C, the length of the compressor coil 

The inductance for a current density entering the compressor coil at lAC will 
be 


Lj = Lj - lAC 


nAC 


( 11 ) 


Then the total inductance L is obtained by adding the different inductances Lj according 
to the rule for parallel circuits. 



( 12 ) 


(13) 


(14) 
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It IS assumed that there is no current flowing into the last element of the compressor 
coil, eliminating the infinite term for i = n This is equivalent to a residual inductance at 
the end of the compressor coil representing the inductance of the backstrap lead to the 
annular electrode Since 


0 < 1 - - < 1 
n 


then 


L 


> 1 


The initial inductance Lj will be larger than the effective inductance The plasma flow 

into the compressor coil will, therefore, reduce the inductance in the equivalent circuit 
from the initial value Lj to the value L 

In the increasing inductance mode, the compressor coil is connected directly to 
the annular electrode Then the inductance in the equivalent circuit will increase as the 
plasma moves into the compressor coil. 

Both configurations were tried The decreasing inductance configuration was 
selected for further investigation after initial expenments indicated that it would be more 
efficient for the acceleration of glass beads with diameters between 100 and 700 microns 
These were attached to a Mylar foil and then placed across the hole in the coil space 
plate 


IV. THEORETICAL INVESTIGATION 


Artisimovitch [8] developed the theory of the rail-type accelerator Other 
investigators have extended this theory to describe the coaxial accelerator The coaxial 
accelerator with a compressor coil can be represented by an equivalent circuit which 
consists of the coaxial accelerator circuit to which is added a mesh with variable plasma 
resistance and a variable inductance simulating the compressor coil An additional switch 
IS incorporated in the mesh to simulate the arrival of the plasma in the coil 


A. The Equivalent Circuit 

In Figure 1 1 the acceleration of the plasma in the coaxial accelerator is 
represented by the variable inductance L 2 Lq and R 2 are the capacitor bank inductance 

and resistance, which are constants The resistance Rj simulates the arrival of the plasma 
at the compressor coil and Lj simulates the compressor coil inductance (see previous 

discussion of compressor coil configuration). The differential equations for this equivalent 
circuit with the additional switch closed are 
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ADDITIONAL SWITCH TO 
SIMULATE ARRIVAL OF 



RESISTANCE 


Figure 1 1 Equivalent circuit for coaxial accelerator with compressor coil 


Uo 


— / (II + 12) dt = Lo — (II + 12) + R2 (II + 12) + R1 II 

Cq -f) dt 


+ ^(Ll II) 
dt 


(15) 


and 


Uq 


1 i d d 

— / (II + 12) dt = Lq - (II + 12) + R2 (II +12) + - (L2 12) 

Cq 0 

(16) 


II is the current in the compressor coil and 12 is the current in the coaxial accelerator. 
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A nondimensional notation is now introduced. With 


COq" 


1 


( 17 ) 


the tune becomes 


t - ^ 

OJQ 


(18) 


and the time derivatives are 


d 

dt 



(19) 


and 


d^ _ 2 


( 20 ) 


The nondimensional currents are 


XI, X2 


11 , 2 

ojQ Uq Cq 


( 21 ) 


the nondimensional inductances are 


X3, X4 


LI, 2 


( 22 ) 
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and the nondimensional resistance terms are 


X5, X6 


Rl, 2 

COQ Lq 


( 23 ) 


Then the differential equations (15) and (16) may be written 


T 

1 - / (XI + X2) dT 
0 


— (XI + X2) + X6 (XI + X2) + X5 XI + — (XI X3) 
dT dT 


(24) 


and 


T 

- f (XI + X2) dT 
0 


d d 

— (XI +X2) + X6(X1 +X2) + — (X2X4) 
dT dT 


(25) 


The initial values at T = 0 are 


XI = 0 


(26) 


and 


X2 = 0 


(27) 


With these values, the first derivatives are determined from equations (24) and (25), 



T=0 


X4 

X3 + X4 + X3 X4 


(28) 
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and 


dX2\ _ X3 

dT / ~ X3 + X4 + X3 X4 

T=0 


(29) 


The differential equations (24) and (25) are then differentiated to eliminate the integrals 


d^ d d d^ 

- (XI + X2) = — - (XI + X2) + X6 — (XI + X2) + — (X5 XI) + — - (X3 XI) 
dT^ dT dT dT2 


(30) 


and 


(XI + X2) = (XI + X2) + X6 4; (XI + X2) + (X2 X4) 

dT'^ dT dT'^ 


(31) 


In these equations X6 is a constant. 


B. Analytical Solution for a Simplified Case 

Equation (31) is subtracted from equation (30) 


^ (XI X3 - X2 X4) + ^ (XI X5) = 0 


(32) 


The result of the first integration is 


d 

— (XI X3 - X2 X4) + XI X5 = C 


(33) 


Using the initial values from equations (26), (27), (28), and (29), 


dXl dX2 

C = X3 - X4 — — = 0 


dT 


dT 


( 34 ) 
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Equation (33) then becomes a first order differential equation for XI if the other 
functions X are known. 


dXl d (In X3) 

— — + XI 

dT dT 



J_ i_ 

X3 dT 


(X2 X4) = 0 


(35) 


The experimental results indicate that during the first half-period of the discharge 
the current through the coaxial accelerator, X2, may be approximated by the undamped 
capacitor bank current with a frequency 


cjQ = , (36) 


and the maximum current at T = 7 t/ 2 becomes 


120 


U 


0 


COq Lq 


or X20 = 1 


(37) 


Then, 


X2 = X20 sin T = sin T 


(38) 


This result imphes that the inductance of the coaxial accelerator does not change, 
X4 = X40 = constant Other approximations are necessary for the functions X3 and X5 
The decrease of the compressor coil inductance is represented arbitrarily as 


X3 


X30 

(1 + T)"’ 


where X30 = X3 at T - 0 


(39) 


and the arnval of the plasma at the compressor coil is approximated by 


X5 


X50 

(1 + T)" 


where X50 - X5 at T - 0 


(40) 
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Then the differential equation (35) for XI becomes 


dXl 

"dT 


+ XI 


d(ln X3) 
dT 


X50 

5 ^ 


n-m 


(1 +T) 


X40 

X30 


(1 + T) cos T 


0 


(41) 


For n = m, the solution for XI is 


XI 


X40 


(1 + T)’^ 



X50 


sin T + exp 



cos T 


(42) 


and if the parameters 


X40 _ L20 

~ ITo 


(43) 


and 


X50 _ RIO 
5^ ~ coQ LIO 


(44) 


are set equal to 1, a very simple expression is obtained. 


XI 


n 

(1 +T) 

2 



cos T + exp (- T)J 


(45) 


The nondimensional functions XI and X2 for the currents are plotted in Figure 12 for 
vanous values of equations (43), (44), and n = 0.5 
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XI AND X2 AS FUNCTIONS OF TIME T 
Figure 1 2. Compressor coil current from equation (44) 


C. Numerical Integration of the Equations for the Equivalent Circuit 

For the numencal integration the differential equations (32) and (33) are 
rewritten 
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and 



+ X2 



[(1 + X3) (1 + X4) - l] 


(47) 


Equations (46) and (47) are integrated in four sequential steps The motion of the 
plasma m the coaxial accelerator is descnbed during the first step The second step begins 
when the plasma enters the compressor coil and ends when the plasma reaches the 
downstream (narrow) end of the compressor coil The next step describes the behavior of 
the system dunng the remaining part of the compression phase which ends at the 
compressor coil current maximum Dunng the fourth step the compressor coil current 
will decrease 

1 Calculation of Step 1 During the first phase of the discharge, when the 
plasma is still mside the coaxial accelerator, there will be no electrical current m the 
compressor coil. To simulate this condition a switch is added to the compressor coil loop 
of the equivalent circuit shown in Figure 1 1 This switch is open as long as the plasma is 
in the coaxial accelerator Step 1 is, therefore, described by the differential equation for 
the coaxial accelerator current which is obtained from equation (33) with XI and 
dXl/dT set equal 0, 


d^X2 

dT^ 




^X6 + 2 



(48) 


The coaxial accelerator inductance X4 changes with the position of the plasma 
dunng the acceleration process. The force acting upon the plasma, which is assumed 

to be a mass m at the position z of the coaxial accelerator, can be determined from the 
partial derivative of the energy W of the configuration with respect to the coordinate m 
which the plasma motion is directed. 


F 


z 



(49) 
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The inductance representing the coaxial accelerator in the equivalent circuit 
changes with the plasma position Z in the coaxial accelerator, 


L2 = L20 + bz 


(50) 


where 


L20 = L2 at T = 0 , z = 0 


(51) 


and b is the gradient of the inductance of the coaxial accelerator 


b 



(52) 


where 


radius of outer electrode 
radius of inner electrode 


(53) 


The force F 2 acting upon the plasma mass m is then 


F -7 = m 


d^Z 
d t^ 


I^b 


(54) 


With the nondimensional coordinate 


X7 = 



(55) 


and the constant 


K7 


b^ Cq^ Uq^ 

2m Lq 


(56) 
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the equation of motion for the plasma is 


d^XV 

dT^ 


k7 (X2)2 


( 57 ) 


The inductance X4 is determined from equation (50) which is divided by Lq for 
nondimensional notation, 


X4 = X40 + X7 


(58) 


with X40 equal to the nondimensional inductance of the coaxial accelerator at T = 0 and 
X7 IS calculated by integrating equations (57) and (58) The initial values for these 
equations at T = 0 are X2 = 0, and 


dX2 _ 1 V 

dT ~ 1 + X4 


(59) 


The plasma motion starts at z = 0, with the velocity 0, therefore, 


^ ,X1 = 0 


(60) 


and 


dX7 

tT 


(61) 


Using these initial values, the differential equations (48) and (57) are integrated until the 
nondimensional coordinate X7 reaches the value attained at the end of the coaxial 
accelerator, X7 = X71. 

2. Calculation of Step 2 At the beginning of Step 2 the switch m the 
compressor coil loop is closed as the plasma moves into the compressor coil In a first 
approximation (see previous discussion of compressor coil configuration), the plasma 
velocity IS set constant and equal to the value at the end of step 1 
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dX7 

dT“ 


( 62 ) 


dX7 \ 

' X7 = X71 


The time needed for the plasma to reach the downstream (narrow) end of the compressor 
coil with the coordinate X72 is then 


T2 = 


XI 2 

dX7 \ 

' X7 = X71 


(63) 


Dunng Step 2, the compressor coil inductance decreases from the initial value 
X31 to a lower value if the configuration operates in the decreasing inductance mode. 
This behavior is approximated by a linear function of the time T, 


X3 


X31 1 


(' 


K1 


T - T1 \ 
T2-T1/ 


(64) 


where 


X31 = X3 at T = T1 


and where a free constant K1 can be used to adjust the theoretical results to fit the 
experiment Another linear function of T and a free constant K2 is used to simulate the 
decreasing resistance between the center electrode and the compressor coil which drops 
from a high initial value X51 to the plasma resistance 


X5 


X51 


(■ 


- K2 


T -T1 \ 
T2 - Tl/ 


(65) 


where 


X51 = X5 at T = Tl 
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It IS assumed that the inductance of the coaxial accelerator does not change after 
the plasma has reached the end of the accelerator at step 1, X4 = X41. The initial values 
for the electrical currents are 


XI =0 and X2 is equal to the. final value of step 1 


( 66 ) 


The first derivatives are determined from equations (24) and (28) at T = T1 Equation 
(28) can be rewritten as 


dXl _ X4 dX2 

dT ~ 3g It" 


(67) 


and equation (24) at T = T1 is 


d d 

1 - / X2dT = — (XI + X2) + X2 X6 + — (XI X3) 

dT dT 


Using equation (67), the first derivative of X2 is obtained 


T=T1 

1 - / X2dT - X2 X6 
dX2 T=0 


dT 

/ 

1 \ 


1 + X4 

1 + — ) 


\ 

. X3/ 


( 68 ) 


(69) 


and the first derivative of XI is 




T=T1 



1 - f X2dT - X2 X6 

dXl 

X4 

T=0 

dF 

X3 



( 70 ) 
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The integral 


T=T1 

f X2 dT (71) 

T=0 


is evaluated dunng Step 1 . Using these initial values, the differential equations (46) and 
(47) are integrated until T = T2 

3 Calculation of Step 3 During Step 3 the inductance X3 and the resistance X5 
remain constant and have the same value as at the end of Step 2. The compressor coil 
current XI and the coaxial accelerator current X2 have initial values which are equal to 
the final values of Step 2 The initial values for the first derivatives of the current 
functions are calculated using equations (24) and (25), which are solved for dXl/dT and 
dX2/dT 


dXl _ X4 [1 - Xll - X6 (XI + X2) - (1 + X4) XI X5] 
"df X3 + X4 + X3 X4 

and 

dX2 _ X3 [ 1 - XI 1 - X6 (XI + X2)] + XI X5 
"dT X3 + X4 + X3 X4 


(72) 


(73) 


T 

where XI 1 = f (XI + X2)dT 
0 


Equations (46) and (47) may now be integrated until the compressor coil current XI has 
reached its maximum. 

4 Calculation of Step 4. Again it is assumed that X3 and X5 do not change with 
time The initial values are obtained in the same way as in Step 3 During this step there 
will be a decompression The integration is terminated after a time corresponding to the 
time recorded on the oscilloscope during the discharge, usually 100 microseconds. The 
computations have been performed at the Leibniz Rechenzentrum of the Bayerische 
Akademie der Wissenschaften at the Technische Universitat Mlinchen, Germany. 
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V. EXPERIMENTAL AND THEORETICAL RESULTS 


The theory developed in the preceding section is compared with expenments 
performed in the Plasma Hypervelocity Range of the Space Sciences Laboratory at the 
George C. Marshall Space Flight Center 


A. Plasma Flow 

The theory can be correlated with a specific experiment through the elements of 
the equivalent circuit shown in Figure 1 1 The inductance of the compressor coil L, and 
the inductance of the capacitor bank with coaxial cables and switch assembly Lq were 
measured pnor to the expenment [11] Throughout the expenments described here, the 
capacitor bank and the switch assembly have not been changed. Then the only unknown 
quantities are the plasma resistance Rj and the constants Kl, K2, and K7 

The plasma resistance has little influence upon the characteristics of the 
equivalent circuit. It may be regarded as part of the impedance Z of the circuit mesh 
with compressor coil inductance Lj , where 


z = yj Rj^ + co^ Li^ 


(74) 


A typical example is 


R1 = 4.8 X 10"^ ohm 


(75) 


and 

LI = 2 X 10'6 hy 

from the experiment described m Table 2 The frequency is 

CO = 9.45 X lO'^ sec‘‘ (76) 

Then the plasma resistance term is three orders of magnitude smaller than the inductance 
term. 
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The increase of the coaxial accelerator inductance Lj in Step 1 is proportional to 
the constant K7 Equation (58) may be rewritten as 


K7 


Uq^ Cq Cq 
2m Lq 


(77) 


The first term of equation (77) is the stored energy per plasma mass m For a coaxial 
accelerator with given geometry and constant plasma mass, K7 changes with the capacitor 
bank charging voltage Uq The results of the computation of the plasma velocity at the 
end of a coaxial accelerator 14 cm long are compared with measured values m Figure 13 
For this example, K7 = 1 0 corresponds to a charging voltage of 18 kV The plasma 
velocity was measured with two light probes located 20 cm apart, one at the end of the 
coaxial accelerator and the other behind the compressor coil The order of magnitude of 
the velocity from theory and experiment is the same The velocity increases with larger 
values of Uq The light probes could only measure the propagation velocity of a luminous 
front, which may not coincide with the major portion of the plasma mass. 

The photographic recording of the oscilloscope trace of the current m the 
compressor coil as well as the capacitor bank voltage indicates the time when the plasma 
left the coaxial accelerator The capacitor bank voltage recordings show irregulanty when 
the compressor coil current begins to rise The arnval time has been plotted in Figure 14 
as a function of the capacitor bank charging voltage Uq For comparison, the recordings 
of the first light detector used for velocity measurement are also shown. The third curve 
in this figure represents the values calculated using the first step outlined in the previous 
section, equations (57) and (48). The expenmental data are listed m Table 3 

1 Capacitor Bank Voltage The capacitor bank voltage measurement was 
described in the section on diagnostics and particle velocity measurement. The voltage 
across the capacitor bank is 


_U 

Uo 


T 

/ (XI + X2) dT 
0 


(78) 


(see equations (24) and (25)) 

The recorded curve is compared with the calculated curve in Figure 15 The 
expenmental curve shows a potential drop approximately 30 microseconds after initiation 
of the discharge This potential drop has also been observed without the compressor coil, 
and it comcides approximately with the transition of the coaxial accelerator current 
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Figure 13 Plasma velocity at the end of the coaxial accelerator. 


through zero. If the coaxial accelerator is still filled with a plasma, then another discharge 
can be initiated, causing a shift in the potential drop There is an indication of the same 
effect after another half-penod of the discharge During these experiments with different 
compressor coil configurations, high-speed motion pictures were tried. A camera with up 
to 20,000 frames per second was used so that the exposure time corresponded roughly to 
one half-penod of the discharge. These motion pictures did not give detailed information 
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Figure 14 Time of arrival of the plasma at the compressor coil. 


about the compressor coil flow, but they showed that there were several subsequent 
accelerations of plasma clouds These observations tend to confirm the explanation given 
for the potential drop observed on the experimental curve The measured frequency of 
the discharge was nearly constant dunng the experiments, varying only by 2 
microseconds The integrations of the equivalent circuit equations with different values 
for the compressor coil inductance showed the same behavior The frequency of the 
discharge is, therefore, mainly determined by the capacitor bank values Cq and Lq • 

2. Compressor Coil Current. The magnitude and the change with time of the 
compressor coil current determine the efficiency of the compressor coil. The experimental 
curve IS compared with the calculated curve in Figure 16 During these experiments, the 
amplitude of the compressor coil current trace was not calibrated. The absolute value of 
the compressor coil current can, therefore, only be estimated from the comparison with 
the theory The theory shows that the magnitude of the compressor coil current depends 
upon the ratio of tlie inductances Lq, L| , and L 2 , all of which are well known It, 
therefore, seems justifiable to use the correlation between theory and experiment to 
estimate the magnitude of the compressor coil current Rewriting equation (30), 


11 = XI coo Uo Co 
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TABLE 3 EXPERIMENTAL AND THEORETICAL DATA 


Data describing the experiment 

Bank capacitance 

Bank inductance 

Bank resistance 

Coaxial accelerator inductance 

Compressor coil inductance 

Coaxial accelerator inductance 
per unit length 

Coaxial accelerator length (for 
plasma acceleration) 

Compressor coil length 

Plasma resistance 

Aluminum foil mass 

Glass particle diameter 

Charging voltage 

Measured particle velocity 

Data used for the theory 


Co 

= 2.7 

X 

10 -^ 

F 

o 

1-1 

= 4.16 

X 

10-’ 

Hy 

R, 

= 79 

X 

10'3 

Ohm 

L, 

= 0.72 

X 

10 '^ 

Hy 

L, 

= 4 43 

X 

10 -’ 

Hy 

b 

= 3.58 

X 

10 -''. 

Hy/m 


20 cm 
13 cm 

R, = 4 8 X 10'3 Ohm 
m = 0.04 g 
D = 0 25 mm 
Uo = 16 kV 
V = 5 3 km/s 


X30 = 

0.4 

K2 

= 0 999 

X40 = 

0 173 

K4 

= 0.324 

X50 = 

1 8 

K5 

= 0.204 

X6 = 

02 

K7 

= 07 

K1 = 

0 5 
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CAPACITOR BANK VOLTAGE AS A FUNCTION OF TIME FOR K7 = 0.7 


Figure 1 5 Capacitor bank voltage 
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COMPRESSOR COIL CURRENT AS FUNCTION OF TIME FOR K7 =0. 7 


Figure 16. Compressor coil current 



With the values given in Table 3, the maximum value of the compressor coil current 
becomes 


11^^^ = 4 08 X 10® XI = 4 96 X 10^ A 

lllclA 


(79) 


The magnetogasdynamic flow m the vicinity of the compressor coil axis can be descnbed 
by the nonviscous magnetogasdynamic momentum equation, 


P 


d v 

IT 


grad p + j X B 


(80) 


The magnetogasdynamic force anses from the interaction of the electric current density 


j = a (E + V X B) 


(81) 


with the magnetic induction B - 



(82) 


The components of the vector product j X B are, m cylindrical coordinates. 


^<p Jz 

J X B =[ Bj. - Jr B^ 

Jr h - H 


(83) 


The radial component will perform the compression of the plasma The current 
density is 


J0 " 


(84) 
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The first term in equation (84) comes from the change of the magnetic field with 


time, 


H = - i V X E 


(85) 


The axial field on the compressor coil axis can be written 


= A(z) I(t) 


( 86 ) 


where A(z) is a function derived from the compressor coil geometry and I(t) is the 
compressor coil current. By integrating equation (85), the azimuthal electrical field from 
the time varying magnetic field becomes 


j M A(z) I(t) 


(87) 


The second term of equation (84) is the product of the axial plasma velocity and 
the radial magnetic field. The latter can be approximated with the assumption that the 
axial magnetic field does not change m the vicinity of the compressor coil axis with the 
radius r. Then 


V H = 0 (88) 


gives an approximation for the radial magnetic field. 


H 


r 


r d A(z) 
2 d z 


I(t) 


(89) 


The azimuthal electrical current density is then 


J0 - - 


o M 
2 


A(z) 


d I(t) 
dt 


+ V, 


d A(z) 
dz 



(90) 
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The first term of the radial force acting upon the plasma is the product of this current 
density with the axial magnetic inductance, 

a [x^ r dl(t) d A(z) 

^ r A(z) I(t) I A(z) I(t) (91) 

As long as Fj. is negative, the force acts towards the axis and the plasma is 

compressed in the coil. The first term is negative during the first quarter period of the 
compressor coil current when the time denvative of the magnetic induction is positive 
The sign of the second term depends upon the direction of the radial magnetic field This 
will be negative approximately the first 80 percent of the coil length The second term m 
the radial component of the magnetogasdynamic force, equation (83), will always be 
negative It is the interaction between the axial current density from the center electrode 
to the compressor coil with the azimuthal magnetic field induced by this current density. 

The azimuthal component of the magnetogasdynamic force produces an azimuthal 
plasma velocity which is small near the compressor coil axis because of the symmetry of 
the configuration. The axial component of the force F^ may not be negative if there is to 

be particle acceleration. The square of the flow velocity in equation (1) indicates that the 
axial term. 


^z Jr J0 ®r ’ 

should always be positive A first approximation would be that these forces are 
transferred without delay into thermodynamic quantities such as pressure and density. The 
compressor coil, therefore, mcreases the density of the plasma flow leaving the coaxial 
accelerator. The expenment shows that the compression is terminated when the 
compressor coil current is at its maximum value. 

The comparison of theory and experiment shows that the equivalent circuit and 
the corresponding theory descnbed in the foregoing section simulate the characteristics of 
the coaxial accelerator with attached compressor coil One of the principal reasons is that 
the behavior of the experiment, as well as the equivalent circuit used for the theory, is 
dommated by the capacitance Cq and the inductance Lq of the capacitor bank, coaxial 
cables, and switch assembly, and the ratio of Lq to the inductances of the coaxial 
accelerator Lj and the compressor coil L, . 

3. Coaxial Accelerator Current The voltage measured at the capacitor bank 
switch assembly can be differentiated to obtain the total current flow, I, +I 2 , to the 
coaxial accelerator and compressor coil There was no separate measurement of the 
coaxial accelerator current dunng the experiments described here because the compressor 
coil was grounded to the coaxial accelerator outer electrode (Fig 3). 
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The theoretical curves for X2 and XI + X2 have been calculated using the theory 
which has been previously shown to agree with the experimental curves of Figures 15 and 
16. They are shown in Figure 17. The third curve shows the current in an equivalent circuit, 
which contains only the capacitor bank elements and the coaxial accelerator The three 
curves are identical during Step 1 of the first and second curves. The third curve is 
obtained by extending the integration of Step 1 until T = 10.0 



CURRENT (XI + X2) AND X2 AS A FUNCTION OF TIME FOR K7 = 0, 7 


Figure 17. Coaxial accelerator current 


B. Particle Acceleration 

The acceleration of the particles placed on a Mylar foil across the narrow end of 
the compressor coil can be used as an indication of the amount of compression of the 
discharge. Figure 18 is a photograph of the beads after they have left the Mylar foil 
They have an initial diameter of 0 8 mm and a velocity of 4 to 5 km/s The velocity 
measurement method has been described previously During the expenments the pressure 
m the range was always kept below 25 microns, and the capacitor bank charging voltage 
and the particle initial diameter were varied. Only sphencal glass beads were used 
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Figure 18. Photograph of accelerated particles in compressor coil plasma flow. 


Figure 19 shows the particle velocity as a function of the capacitor bank charging 
voltage for three initial particle diameters. The dashed curve in these figures is obtained 
from the assumption that after the acceleration, kinetic energy of the particles with the 
mass mp is proportional to the energy stored in the capacitor bank, CqUo^/ 2, multiplied 

by the particle cross section S and a coefficient C, representing the energy transfer 
efficiency 


m. 


,2 = 


= c, s 


n 2 IT 2 
'-'0 ^0 


(93) 


Then the velocity of a particle with the initial diameter D is given by 


V = Cj 



where C 2 


y 3 c. cq 

“ ^GLASS 


(94) 
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Figure 19 Particle velocity as function 
of capacitor bank charging voltage. 


In the case of the 700-micron 
diameter particle, the straight Ime 
representing this theory could easily be 
drawn through the experimental results The 
theoretical line for another particle initial 
diameter will have an inclination which is 
related to the 700-micron curve inclination 
by 


tan a. 


tan Of^QQ 



(95) 


The inclinations of the curves in the 
figures for the 250-micron and the 
125-micron particles were calculated from 
this formula The absolute level of these 
curves was adjusted to the 700-micron 
experiments The constant m equation 
(94) was determined from the line 
approximating the 700-micron results The 
lines for 250 microns and 125 microns can 
then be computed with equation (94) In 
Figure 20 the final velocity is plotted over 
the particle initial diameter The theoretical 
curves are also calculated from equation 
(94) 


The difference between the velocities measured in different experiments with the 
same initial particle size mcreases with large capacitor bank charging voltage, very 
probably due to particle ablation dunng the interaction with the plasma 


VI. DESCRIPTION OF THE COMPUTER PROGRAM* 


The differential equations are integrated numerically in four sequential steps by 
the SUBROUTINE HPCG The equations for each step are given m the SUBROUTINES 
FCTl, FCT2, FCT3, and FCT4, which are called by SUBROUTINE HPCG. The output is 
performed through the SUBROUTINE OUTP 


1. A listing of the computer program can be found in the Appendix. 
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PARTICLE VELOCITY [10’m/s] 



Figure 20 Particle velocity as function of the particle initial diameter 










In the program some names are used which are different from those used in the 
Theoretical Investigation section to avoid changing most of the integration statements in 
SUBROUTINE HPCG, which was provided by J. Steuerwald^ . The corresponding names 
were 


Independent variable T ->• X 

Dependent variable XI, X2, . XI 1 ^ Y(l), . . Y(1 1) 

Constants Kl, K2, .. K7^VK(1), .. VK(7) 


A. Input Data 

The input data are given by the MAIN PROGRAM The constants VK(i) and 
VK(2) were introduced and described in equations (64) and (65), VK(7) in (56), and 
VK(4) and VK(5) represent the coordinates of the ends of the coaxial accelerator and the 
compressor coil, respectively. The mitial values Y(l) . Y(ll) for the computation of 
starting values for the first integration-step are Y(l), Y(2), Y(7), Y(8), Y(10) and Y(ll) 
are equal to zero at T = 0, Y(3), Y(4) and Y(5), Y(6) represent the nondimensional 
values for the inductances LI, L2 and the resistances Rl, R2. Y(9) is obtained from 
equation (59) After the output statements for the listing of the input data, several 
parameters of the integration such as the bounds of the integration interval and the error 
bound and the initial increment of the independent variable X are furnished to the 
integration program. 


B. Integration 

It was not necessary to use different error weights The initial value Y(3) for the 
coil inductance is stored on “core” for some calculations during STEP 3. If the plotter is 
on-line, an index (IND) is used to count the calculated elements which are moved to the 
plotter-storage-array . 

With these specifications, SUBROUTINE HPCG begins the integration of the first 
equation system which is defined in SUBROUTINE FCTl. DERY(l) . DERY(ll) are 
the first denvates of the independent variables DERY(9) is obtained by equations (28) 
and (29), DERY(IO) by (57), and DERY(ll) is described in (73) The parameter KKX 
indicates the number of the integration step, the steps are 1, 2, 3 or 4 and are used by 
SUBROUTINE OUTP in a control statement for termination HPCG integrates the 
equations until there is satisfactory accuracy. The calculated values are then transmitted 
to SUBROUTINE OUTP, where some of them are changed from nondimensional to real 
values 


2. J Steuerwald, Instituet fur R. Plasmaphysik, Garching, Germany, private communication 
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Real Time 


RT = X/WO = X 1 .06 (10'5 sec) 


Real Place RP = Y(7) Lo/b = Y(7) 1 16 (m) 

Real Velocity RV = Y(10) Lo WO/b = Y(10) 1 1 (10^ m/sec) 


The electncal currents Y(l), Y(2), Y(l)+Y(2) and the voltage l-Y(ll) are listed as 
nondimensional values. Time X, voltage U, and the electncal currents Y(l) and Y(2) are 
moved to a plotter-storage-array which is called at the end of the MAIN PROGRAM. 

After each integration, a control statement is called which is specified by the 
parameter KKX STEP 1 is terminated by means of one of these control statements when 
the plasma has reached the end of the coaxial accelerator, that means RP = VK(5). 

The mitial values for STEP 2 are defined at the end of STEP 1. V2 represents 
equation (49), Y(8) and Y(9), the current derivatives, are determined by equations (55) 
and (56) m SUBROUTINE FCT2, DERY(3) and DERY(5) are the functions for the 
decreasing inductance and resistance introduced in equations (64) and (65). DERY(8) and 
DERY(9) are defined by equations (30), (31), (46) and (47). The integration of STEP 2 
is continued until the plasma reaches the end of the compressor coil; RP is then equal to 
or greater than VK(4) 

After termination of STEP 2, SUBROUTINE OUTP initializes the integration of 
the equation system given in SUBROUTINE FCT3 The initial values for Y(8) and Y(9) 
are obtained from equation (58) DERY(8) and DERY(9) in SUBROUTINE FCT3 
correspond to equations (30) and (31), with DERY(3) = DERY(4) = DERY(5) = 0. This 
system is mtegrated as long as the first derivative Y(8) of the compressor coil current is 
negative. The integration of STEP 3 is terminated by the change of the sign of Y(8) 

The integration of STEP 4 begins with the initial values for Y(8) and Y(9), which 
are obtamed by equations (28) and (29). The system of equations in SUBROUTINE FCT4 
IS the same as in SUBROUTINE FCT3. The integration is continued until the 
independent variable X has reached the upper bound of the integration interval which 
was specified in the first section of the MAIN PROGRAM 


C. Example of Calculations 

The computer program described m this section may be used to investigate the 
influence of different parameters on the equivalent circuit representing the coaxial 
accelerator and attached compressor coil One of the parameters which can be varied is 
the coefficient K7 which is defined m equation (56). The parameter K7 represents the 
energy stored in the capacitor bank In Figure 21 the curves for the nondimensional 
capacitor bank voltage U/Uq, the compressor coil current XI, and the coaxial accelerator 
current X2 are shown. The corresponding values of the variables at T = 0 and the 
parameters are listed in Table 3, and K7 = VK7 is vaned as follows 
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Figure 2 1 Voltage and current curves for different values of the constant K7 



VK7 = 0.543 = 14 KV 
VK7 = 07 = 16 KV 
VK7 = 0.886 = 18KV 
VK7 = 1.09 = 20 KV 


The value of VK7 was varied since VK7 ~ Uq, as shown in equation (77). Figure 21 
shows that the compressor coil current XI decreases as K7 increases if all other 
parameters are held constant. 

Finally, Figure 22 shows the variation of XI, X2 and U with the plasma 
resistance X5. The current XI increases if all other parameters have the value shown in 
Table 3. 


VII. CONCLUSION 


The expenmental results show that a coaxial accelerator with a compressor coil 
produces a higli density plasma, which can be used to accelerate glass particles for 
micrometeoroid simulation The discharge of energy stored in the capacitor bank can be 
investigated theoretically with an equivalent circuit. The theoretical and experimental 
curves for the electrical current flowing in the coaxial accelerator and in the compressor 
coil agree withm the accuracy of the theoretical approach and the measured values for 
the equivalent circuit components The photographs of the compressor coil flow and the 
particle acceleration, together with the high velocities obtained, indicate that the 
compressor coil confines the plasma dunng the time when the compressor coil current is 
increasing 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, July 1974 
502-21-28-0000 
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Figure 22. Voltage and current curves for different values of the plasma resistance X5 
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COMPUTER PROGRAM LISTING 
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PROGRAM HYPERV 


C 
C 
C 

c integration by A steps 

C 

COMMON KKX,V1,V2,V3.V(.,V5,U4,W1,W2,W5 
COMMON VK(8) 

common IND,TPLOT(300),UPLOT(300),Y1PLOT(300),Y2PLOT(300) 
dimension Y (1 1) ,DERY h 1) .PRMT (5) ,AUX (16, 1 1) 

c 

C DATA-INPUT 

c 

c constants 

Vk (1) =0.5 
VK (2) =0.999 
Vk U) =0.27 
VK (5) =0.14 
VK (7) =. . . . 

INITIAL VALUES (STEPI) 

X = 0.0 
Y<1)=0.0 
Y(2)=0.0 
C Y(3)=.... 

C Y (4) =. . . . 

C Y (5) = . . . . 

C 

READ (5,1) VK(7) ,Y(3) ,Y(4) ,Y(5) 

1 FORMAT (4(F10.4)) 

Y (6) =0.2 
Y (7) =0.0 
Y (8) =0.0 

Y(9)=1./(1.*2.*Y(41) 

Y (10) =0.0 
Y(11)=0.0 
C 

C OUTPUT-FORMATS 
C 

WRITE (6,4) 

WRITE (6,5) (Y(I),I=1,7) 

WRITE (6,6) 

WRITE (6,7) (Y (I) , 1=8, 1 1) 

WRITE (6,8) 

WRITE (6,9) VK (1 ) , VK (2)- , VK (4) , VK (5) , VK (7) 

WRITE (6,11) 

WRITE (6,12) 

WRITE (6,17) 

C 

4 FORMAT (1H1///,2X,29HINITIAL VALUES AND C 0 N S T A N T S : , / , 2 X , 29H 

1 ,///,5X,4HY(1),8X,4HY(2),8X,4HY(3),8X,4HY(4), 

28X,4HY(5) ,8X,4HY(6) ,8X,4HY(7)) 

5 FORMAT (7(2X,E10.3) ) 

6 FORMAT (/,5X ,4HY (8) ,8X,4HY (9) ,8X,5HY (10) ,7X ,5HY (1 1) ) 

7 FORMAT (4(2X,E10.3) ) 

8 format (/,5X,5HVK(1),7X,5HVk(2),7X,5HVk(4),7X,5HVk(5),7X,5HVK(7)) 

9 FORMAT (5(2X,E10.3)) 

11 FORMAT (////, 2X. 6HSTEP1 ;) 

12 FORMAT (2X,6H ) 

17 F0RmAT(//,4X,5HT (R).8X,4HY(1),8X,4hY(2),6X,9HY(1)*Y(2),5X,4HY(7), 
18X,5HY(10) ,6X,7H1-Y(11) ,/) 


51 



c 

c 

C THE DIFFERENTIAL EQUATIONS ARE INTEGRATED BY SUBROUTINE HPCG 

C 

C NUMBER OF DIFF. EQUATIONS 

NDIM=1 1 
C 

C LOWER BOUND OF THE INTERVAL 

PRMT (1) »0. 

C 

C UPPER BOUND OF THE INTERVAL 

PRMT (2) =10. 

C 

c initial increment of the independent variable X 
PRMT (3) «0.05 
C 

C UPPER ERROR BOUND 

PRMT (A) =1 .0E-03 
C 

C SUBROUTINES CALLED BY SUBROUTINE HPCG 

EXTERNAL OUTP 
EXTERNAL FCT1 
EXTERNAL FCT2 
EXTERNAL FCT3 
EXTERNAL FCTA 
C 

CORE = Y (3) 

C 

C ERROR weights 

DO 20 1=1,11 

20 DERY (I) =1 .0/1 1 .0 
C 

C PLOTTER-INDEX! IND 

IND = 1 
C 

C INTEGRATION STEP1 

C 

call HPCG (PRMT , Y , DERY.NDIM, IHLF. FCT1 , OUTP, AUX) 

C END INTEGRATION STEPI 

C 

WRITE (6,13) 

13 FORMAT (//,2X,6HSTeP2:) 

WRITE (6,12) 

WRITE (6,17) 

C 

C INITIAL VALUES (STEP2) ARE DEFINED AT THE END OF STEPI 

X = PRMT (1) 

V1=X 

V2= (VK (A) -VK (5) ) /Y (10) 

V3 = Y (3) 

VA = V (A) 

V5=Y (5) 

Y(8)=Y(A)*(1.-Y(11)-Y(6)*Y(2))/(Y(3)+Y(A)+Y(3)*Y(A)) 

Y(9)=Y(3)*(1.-Y(1D-Y(6)*Y(2))/(Y(3)+Y(A)*Y(3)*Y(A)) 

C 

C ERROR WEIGHTS 

DO 21 1=1,11 

21 DERY(I)=1. 0/11.0 
C 

C INTEGRATION STEP2 
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c 

call HPCG (PRMT , Y, DERYiNDIM, I HLF, FCT2.0UTP, AUX) 

C END INTEGRATION STEP2 

C 

WRITE (6,U) 

U FORMAT (//.2X.6HSTEP3:) 

WRITE (6,12) 

WRITE (6,17) 

C 

C INITIAL values (STEP3) ARE DEFINED AT THE END OF STEP2 

X=PRMT(1) 

Y(8)'>(Y(A)*(1.0-Y(l1)-Y(6)*(Y(1)+Y(2)))"(1.0 + Y(i.))*Y(5)*Y(l))/(Y(3 
1 ) +Y (A) +Y (3) *Y (A) ) 

Y(9)»(V(3)*(1.0-Y(i1)-Y(6)*(Y(1)*Y(2)))*Y(5)*Y(1))/(y(3)+Y(A)*Y(3) 
1*Y (A) ) 

C 

C ERROR WEIGHTS 

DO 22 1=1,11 

22 DERY (I) =1 .0/1 1 .0 
C 

C INTEGRATION STEP3 

C 

CALL HPCG (PRMT ,Y, DERY, NDIM, IHLF, FCT3,0UTP, AUX) 

C END integration STEP3 

C 

WRITE (6,15) 

15 FORMAT (//,2X,6HSTEPA:) 

WRITE (6,12) 

WRITE (6,17) 

C 

C INITIAL values (STEPA) ARE DEFINED AT THE END OF STEP3 

X°PRMT (1) 

Y (3) =CORE 

Y(8)=(V(A)*(1.0-Y(i1)-Y(6)*(Y(1)*Y(2)))-(1.0+Y(a))*Y(5)*Y(1))/(Y(3 
1) *Y (A) *Y (3) *Y (A) ) 

Y(9) = (V(3)*(1.0-Y(ii)-Y(6)*(Y(1)+Y(2)))+Y(5)*YM))/(Y(3)+Y(A)*Y(3) 
1*Y(A)) 

C 

C ERROR WEIGHTS 

DO 23 1=1,11 

23 DERY (I) =1.0/11.0 
C 

C INTEGRATION STEPA 

C 

call HPCG (PRMT. Y. DERY, NDIM, IHLF.FCTA.OUTP, AUX) 

C END INTEGRATION STEPA 

C 

WRITE (6,16) 

16 FORMAT (////) 

C 

C plot-program begin (BENSON-PLOTTER) 

c 

c 

C PLOTTER-SCALE-FACTORS 

TPLOT (IND)=0.0 
TPl0T(IND*1)“1.0 
UPL0T(IND)»-1.0 
UPLOT ( I ND+1 ) =0. 25 
Y1 PLOT ( IND) »-1 .0 
Y1PL0T(IND+1)»0,25 
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YZplOT ( 1 ND) =-1 .0 
Y2P LOT ( I ND+1) =0.25 
C 

CALL PLOT (10. ,7. , -3) 
call FACTOR(2.0) 

call ax 1 S (0. , a. . 6HT IME T . -6 . 1 0 . , 0 . , T P LO T ( I N D ) , T P LO T ( I N D+ 1 ) ) 

CALL AX I S (0. , 0. , 9HV0LTAGE U , 9 . 8 . , 9 0 . , U P LO T (I N D ) i U P L 0 T ( I ND+ 1 ) ) 

CALL AX 1S(-1.,0..12HEL. CURRENT Y,12,8.,90.,Y1PLOT(IND) ,Y1PL0T(IND+ 
ID) 

call LINE(TPlOT(1).UPlOT(1),IND-1, 1,10,1) 
call LINE(TPL0T(D .Y1PL0T(i),1ND-1,1,0,1) 
call L1NE(TPL0T(1).Y2PL0T(1),1ND-1,1,10,2) 

c 

call SYMBOL(0.2,2.0.0.21»11HHYPERV-HPCG,0.0,11) 

CALL SYMBOL (0.2, 1 .5,0.14, SHU ,0.0,5) 
call symbol (999., 999. ,0.08, 1 ,0. 0,-1) 

CALL SYMBOL (0.2, 1 .2, 0.14, SHY (1) ,0.0,5) 

CALL SYMB0L(999. , 999. ,0. 14, 1H-, 0.0,1) 

CALL SYMBO L (0 .2, 0.9 ,0. 1 4,5HY (2) ,0.0,5) 

call symbol ( 999 . . 999 , , 0 . 08 , 2 . 0 . 0 ,- 1 ) 

C 

CALL PLOT (14. ,0. , -3) 
call FACTOR(I.O) 
call Plot (0.0, -7. 0.-3) 
call PL0T(0.,0.,4) 

STOP 

c 

c plot-program end 
c 

END 

c 

SUBROUTINE HPCG ( PrMT , Y , D ER Y , ND 1 M , 1 H L F , FC T , 0 U T P , AU X ) 

c - - 

c 

SUBROUTINE OUTP ( X . Y , DE R Y , I H L F , ND I M , PRMT ) 

C - 

c 

COMMON KKX,V1,V2,V3.V4,V5,U4,Wl,W2,W5 
COMMON VK(8) 

COMMON 1ND,TPlOT(300),UPlOT(300),YiPlOT(300),Y2PLOT(300) 
dimension Y (1 1) , DERY (1 1) , PRMT (5) 

C 

SS = Y (1 ) +Y (2) 

U = 1 .0-Y(11) 

RT = X*1 .06 
RP = Y (7) *1.16 
RV=Y (10) *1 . 1 
PRMT (1) =X 
C 

UPLOT (IND) »U 
TPLOT (IND) =RT 
Y1PlOT(IND)"Y(1) 

Y2PL0T ( IND) =Y (2) 

C 

WRITE (6,30) RT,Y(1) ,Y(2) ,SS,RP,RV,U 
30 FORMAT (7 (2X.E10.4) ) 

IND=IND+1 

C 

GO TO (100,200,300.400) ,KKX 
C 
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1 00 

IF (RP . 
RETURN 

GE. 

VK(5) ) 

GO 

TO 

500 

200 

IF (RP . 
RETURN 

GE. 

VK(4) ) 

GO 

TO 

500 

300 

IF (Y(8) 
RETURN 

.GE 

. 0.) 

GO 

TO 

500 

400 

RETURN 






500 

PRMT (5) = 

1 . 






RETURN 

END 

C 

SUBROUTINE FCT1 (X.YiDERY) 

C 

c 

COMMON KKX,V1.V2,V3.V4,V5,U4,W1,W2,W5 
COMMON VK(8) 

DIMENSION Y (1 1) ,DERY (1 1) 

C 

DERY(I) =0. 

DERY(2)=Y(9) 

DERY (3) =0. 

DERY(4)=»Y(10) 

DERY (5) *0. 

DERY (6) =0 . 

DERY(7)=Y(10) 

DERY (8) =0. 

DERY(9)=(-Y(2)*(VK(7)*Y(2)*Y(2)+1.)-Y(9)*(Y(6)+2*Y(10)))/(1.+Y(4)> 
DERY (10) «VK (7) *Y (2) *Y (2) 

DERYM1)=Y(2) 

C 

KKX>=1 

C 

RETURN 

END 

C 

SUBROUTINE FCT2 (X.Y.DERY) 

C 

C 

COMMON KKX,V1,V2,V3.V4,V5.U4,W1,W2,W5 
COMMON VK(8) 

DIMENSION Y(n) .DERYdl) 

C 

DERY (1) -Y (8) 

DERY (2) =Y (9) 

DERY (5) =-Vk (1) *V3/V2 
DERY (4) =0, 

DERY (5) '-V5*VK (2) /v2 
DERY (6) =0. 

DERY (7) =Y (10) 

DERY(8)=(-(1.0+Y(4))*(Y(8)*(Y(6)*Y(5)-2*V3*VK(1)/V2)+Y(9)*Y(6)+Y(1 

1) *(1.0-V5*VK(2)/V2)4.y(2))+Y(6)*Y{8)+V(6)*Y(9)+Y(1)+Y(2))/((1.0+Y(3 

2) ) * (1 ,0*Y (4) ) -1 .0) 

DERY(9)=(-(1.0+Y(3))*(Y(6)*Y(8)*Y(6)*Y(9)+Y(l)+Y(2))+Y(8)*(Y(5)+Y( 

l6)-2*V3*VK(1)/V2)*Y(«)*Y(9)+Y(1)*(1.0-y5*VK{2)/V2)+Y(2))/((1.0*Y(3 

2))*(1.0+Y(4))-1.0) 

DERY (10) =0. 

DERY (11) =Y (1) *Y (Z) 

C 

KKX = 2 
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c 

RETURN 

END 

C 

SUBROUTINE FCT3 (X.Y.DERY) 

C 

c 

COMMON KKX,V1.V2,V3,V4.V5.U4,W1.W2,W5 
COMMON VK(8) 

DIMENSION YdD.DERYdl) 

C 

DERY (1) =Y (8) 

DER Y (2) =Y (9) 

DERY (3) =0. 

DERY (4) =0. 

DERY (5) =0. 

DERY (6) =0. 

DERY (7) =Y d 0) 

DERY{8)=«(-d.0 + Y(4))*(Y(8)*(Y(6)+Y(5))+Y(9)*Y(6)+Yd)+Y(2))+Y(6)*Y 
1(8)+Y(6)*Y(9)+Y(1)+Y(2))/(d.O + Y(3))*d.O + Y(4))-1.0) 

DERY(9) = (-d.0 + Y(3>)*(Y(6)*Y{8)*Y(6)*Y(9)+Yd)+Y(2))+Y(8)*(Y(5)+Y( 
16) ) + Y (6) *Y (9) + Y d) + Y (2) ) / ( d .0 + Y (3) ) *d .0 + Y (4) ) -1 .0) 

DERY (10) -0. 

DERY (1 1) =Y (1) +Y (2) 

C 

KKX = 3 
C 

RETURN 

END 

C 

SUBROUTINE FCT4 (X.Y.DERY) 

C 

c 

COMMON KKX,V1.V2,V3,V4,V5.U4,M1,W2,W5 
COMMON VK(8) 

dimension Y (11) .DERY (11) 

C 

DERY (1) »Y (8) 

DERY (2) =Y (9) 

DERY (3) =0. 

DERY (4) =0. 

DERY (5) =0. 

DERY (6) «0. 

DERY(7) =Y(10) 

DERY(8)=(-(1.0+Y(41)*(Y(8)*(Y(6)*Y(5))+Y(9)*Y(6)+Y(1)+Y(2))+Y(6)*Y 
1 (8) *Y (6)*Y (9) +Y (1) +Y (2) ) / ( (1 .0 + Y (3) )*(1 .0 + Y (4) ) -1 .0) 

DERY(9)*(-(1.0 + Y(3))*(Y(6)*Y(8)*Y(6)*Y(9)+Yd)+Y(2))*Y(8)*(Y(5)+Y( 
16) ) + Y (6) *Y (9) + Y (1) + Y (2) ) / ( d .0*Y (3) ) *d .O + Y (4) ) -1 .0) 

DERY (10) =0. 

DERY (1 1) »Y (1) +Y (2) 

C 

KKX = 4 

RETURN 

END 
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